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Abstract 

We consider non-self adjoint perturbations of a self-adjoint /i-pseudodiffer- 
ential operator in dimension 2. In the present work we treat the case when 
the classical flow of the unperturbed part is periodic and the strength e of 
the perturbation satisfies h s ° < e < e for some 5q e]0, l/2[ and a sufficiently 
small eo > 0. We get a complete asymptotic description of all eigenvalues in 
certain rectangles [—1/(7,1/(7] + ie[Fo — 1/(7, Fq + 1/(7]. In particular we are 
able to treat the case when e > is small but independent of h. 12 



1 Introduction 

This paper is a continuation of [14], where A. Melin and the author observed that 
for a wide and stable class of non-self adjoint operators in dimension 2 and in the 
semi-classical limit (h — > 0), it is possible to describe all eigenvalues individually 
in an /i-independent domain in C, by means of a Bohr-Sommerfeld quantization 
condition. Notice that the corresponding conclusion in the selfadjoint case seems to 
be possible only in dimension 1, or in higher dimensions under strong (and unstable) 
assumptions of complete integrability. In [14] we exploited the absence of small 



*UMR 7640, CNRS 

keywords: Eigenvalue, non-self adjoint 

2 MSC 2000: 32A25, 34M99, 35P20, 35Q40, 37G99 



1 



denominators to get a geometric analogue of the KAM-theorem via methods of non- 
linear Cauchy-Riemann equations and got a corresponding result at the level of 
operators 

In the present work we make another step by studying small perturbations, 
roughly of the form P + ieQ, of a selfadjoint /i-pseudodifferential operator P whose 
associated classical flow is periodic. We will here be particularly interested in the 
case of a small but fixed e, but our methods allow us to let e vary in an interval 
[h 5 °, e ] where e > is sufficiently small and <5 G]0, 1/2 [ is arbitrary. 

From the point of view of applications, it is clear that even smaller perturbations 
are of a considerable interest and as another step, Hitrik and the author [9] studied 
the same problem as in the present paper, but in the parameter range h <C e < h 5 for 
every fixed 5 > 0. When the subprincipal symbol vanishes we could even treat the 
range ti 2 <C e < h s . Actually with M. Hitrik, we plan a whole series of works devoted 
to small perturbations of non-selfadjoint operators in two dimensions. Among other 
things we plan to treat the case when the classical flow of the unperturbed operator 
admits certain invariant Lagrangian torii with a diophantine condition. (Another 
work ([16]) deals with resonances generated by a closed hyperbolic trajectory and 
can be viewed as descendant of the pioneering work of M. Ikawa [10] about scattering 
poles for two strictly convex obstacles.) 

The methods in [9] are partly more traditional and rely on reduction by averaging 
to a one-dimensional problem in the spirit of [21, 5, 4, 6, 11]. Such a reduction does 
not seem possible here and the problem remains 2-dimensional. In general, we 
have been motivated by recent progress around the damped wave-equation ([13], 
[2], [17], [8]), as well as the problem of barrier top resonances for the semi-classical 
Schrodinger operator ([12]) where more complete results than the corresponding 
ones for eigenvalues of potential wells ([18], [3], [15]) seem possible. Eventually we 
also hope to apply our results (though not specifically the ones of the present work) 
to the distribution of resonances for a strictly convex obstacle in R 3 . See [20] and 
references given there. In the case of analytic obstacles, much more can probably 
be said, especially in dimension 3 (and 2). 

The present work was undertaken before the start of [9], but the latter work is 
now completed, so we can take advantage of many of the arguments there, even 
though the main step here will be quite different. 

Let M denote R 2 or a compact real- analytic manifold of dimension 2. 

When M = R 2 , let 

P e = P(x,hD x ,e;h) (1.1) 

be the /i-Weyl quantization on R 2 of a symbol P(x,l;,e; h) depending smoothly on 
e G neigh (0, R) with values in the space of holomorphic functions of (x, £) in a 
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tubular neighborhood of R 4 in C 4 , with 

\P(x,S,e;h)\<Cm(Re(x,S)) (1.2) 

there. Here m is assumed to be an order function on R 4 , in the sense that m > 
and 

m(X) <C (X- Y) No m(Y), X, Y e R 4 . (1.3) 

We also assume that 

m>l. (1.4) 

We further assume that 

oo 

P(x,Z,e;h)~J2Pj,e&0h j , h^O, (1.5) 

3=0 

in the space of such functions. We make the ellipticity assumption 

\po,e&0\ > ^m(Re(x,Z)), \(x,£)\ > C, (1.6) 

for some C > 0. 

When M is a compact manifold, we let 

Pe= o, a , e {x;h){hD x ) a , (1.7) 

|a|<m 

be a differential operator on M, such that for every choice of analytic local coordi- 
nates, centered at some point of M, a Q)e (x; h) is a smooth function of e with values 
in the space of bounded holomorphic functions in a complex neighborhood of x — 0. 
We further assume that 

oo 

a a , e (x; h) a a ,e,j(x)h 3 , /i -> 0, (1.8) 
i=o 

in the space of such functions. The semi-classical principal symbol in this case is 
given by 

P0,e(z, = Yl a «M X )C (1-9) 

and we make the ellipticity assumption 

|po(z,OI > ^<0 m , foO e T*M, |e| > C, (1.10) 
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for some large C > 0. (Here we assume that M has been equipped with some 
Riemannian metric, so that |£| and (£) = (1 + l^l 2 ) 1 ^ 2 are well-defined.) 
Sometimes, we write p e for p 0<e and simply p for p . Assume 

P e=0 is formally selfadjoint. (1.1 1) 

In the case when M is compact, we let the underlying Hilbert space be L 2 (M, fi(dx)) 
for some positive real-analytic density /i(dx) on M. 

Under these assumptions P e will have discrete spectrum in some fixed neighbor- 
hood of G C, when h > 0, e > are sufficiently small, and the spectrum in this 
region, will be contained in a band \lmz\ < 0(e). The purpose of this work as well 
as of [9] and later ones in this series, is to give detailed asymptotic results about the 
distribution of individual eigenvalues inside such a band. 

Assume for simplicity that (with p = p t =o) 

T := p _1 (0) D T*M is connected. (1.12) 

Let H p = p'^ ■ ^ — p' x ■ J| be the Hamilton field of p. In this work, we will always 
assume that for E e neigh (0, R): 

The H p -Row is periodic on T E := p~ l (E) n T*M with (1.13) 
a period T(E) > depending analytically on E. 

(In Section 2 we recall how this assumption follows from a seemingly weaker one.) 
Let q = \ (§-J t=Q Pe, so that 

p e =p + ieq + 0(e 2 m), (1.14) 
in the case M = R 2 and p t = p + ieq + 0(e 2 (^) m ) in the manifold case. Let 

I rT(E)/2 

(q) = / q o exp (tH p ) dt on p~\E) fl T*M. (1.15) 

1 \rj) J-T(E)/2 

Notice that p, (q) are in involution; = H p (q) =: {p, (q)}. As in [9], we shall see 
how to reduce ourselves to the case when 

p € =p + te(q) + 0(e 2 ), (1.16) 

near p~ 1 (0) fl T*M. An easy consequence of this is that the spectrum of P e in 
{z e C; |Rez| < 5} is confined to ] - 5,5[+ie\(Req) minfi - o(l), (Req) maX:0 + o(l)[, 
when S,e,h^ 0, where (Reg) minj0 = min p -i (0)nT » M (Re g) and similarly for (g) maXi0 - 
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We will mainly think about the case when (q) is real-valued but we will work under 
the more general asumption that 

Im (q) is an analytic function of p and Re (q), (1-17) 

in the region of T*M, where \p\ < 1/(9(1). 

Let F G [(Reg) m i nj o, (Reg) ma x,o]- The purpose of the present work is to deter- 
mine all eigenvalues in a rectangle 

] " WV W) [+ie]Fo ~WY Fo + W) 1 (L18) 

for 

h So < e < e , (1.19) 

for h 5 ° < e < e with < 5 < 1/2 and e sufficiently small but fixed. We assume 
that 

T(0) is the minimal period of every iJ p -trajectory in A ,f , (1-20) 

where 

Ao,F := {p e T*M; p{p) = 0, Re (q)(p) = F }, (1.21) 
We also assume that 

dp, dRe{q) are linearly independent at every point of A 0i f - (1-22) 

This implies that every connected component of Ao,f is a two-dimensional La- 
grangian torus. For simplicity, we shall assume that there is only one such compo- 
nent. Notice that in view of (1.20), the space of closed orbits in p _1 (0) fl T*M; 

E : = (p-\0) f]T*M)/ ~, 

where p ~ /i if p — exptH p p for some t G R, becomes a 2-dimensional symplectic 
manifold near the image of A 0i f , and (1.22) simply means that Re (q), viewed as 
a function on E, has non-vanishing differential. The image of A 0i f is just a closed 
curve. 

In [9] (for e in the range h <^ e < h 5 and sometimes /i 2 <C e < h s , W5 > 0) we 
also studied the case when F is a non-degenerate extreme valule of (q) on E. It 
would be interesting to see to what extent that can be done for e in the range (1.19). 

As in [14], the analyticity assumptions seem to be quite essential at least in 
the case of fixed e. Indeed one is naturally led to work in modified Hilbert spaces 
defined by introducing microlocal exponential weights in the spirit of [19, 7, 14, 9], 
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and there are closely related Fourier integral operators with complex phase some of 
which have associated complex canonical transformations that are e-perturbations 
of the identity. 

The plan of the paper is the following: 

In section 2, we make the geometrical work and construct invariant torii close to the 
real domain. This allows us to construct a complex canonical transformation which 
reduces p to a function on the cotangent space on the standard 2-torus, which is 
independent of the space-variables. 

In section 3 we perform further reductions for the whole operator and obtain a 
complete asymptotic description of all the eigenvalues of P e in a rectange of the 
form (1.18). This is still somewhat formal, but 

in section 4, we introduce a global Grushin problem, and verify that the formal eigen- 
values of the preceding section coincide modulo 0(h°°) with the actual eigenvalue 
in a rectangle (1.18). 

2 Geometric reductions 

We use the notation and general set-up of the introduction. Thus let p denote 
the semi-classical principal symbol of the unperturbed operator. As a warm-up we 
recall how the assumption (1.13) follows from a seemingly weaker assumption. Thus 
replace (1.13) by the assumption that for some a > 0, every point p G — a, «[) 
belongs to a closed ifp-trajectory •y(p) with period T(p) > 0, depending continously 
on p. Also assume dp ^ on IV Then, 

1) If 7(p) G T E is the T(p)-periodic if p -trajectory passing through p G p~ l (E), then 
the action I{j(p)) = f,-. £dx only depends on E but not on p. 

2) We have the same conclusion for the period T(p) and hence (1.13) holds. 

Indeed, consider first two trajectories 7(po),7(/°i) an d take an intermediate fam- 
ily j(p s ), < s < 1, depending continuously on s, so that the union of the j(p s ) is 
a two-dimensional manifold T C p" 1 (E). Notice that a\ r = 0, since H p is tangent 
to T and belongs to the radical of the restriction of a to p^ 1 (E). Hence by Stokes' 
formula, 




This shows 1). As for 2), let 7^ C T^; be a smooth family of if p -periodic curves 
with period = T(je)- Let T = UeoKekeiJe and let v be a vector field on T, with 
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u(p) = 1. Let t be a multivalued time variable on T, so that i/ p t = 1. Then we 
claim that 

cr| = dp A dt = d(pdt) : 
On the one hand, (a, v A H p ) = (dp, v) — 1 and on the other hand 

<*A*,,,Aiw = d*(<*$ = 

since the diagonal elements of the matrix are equal to 1, and the claim follows. 
By Stokes' formula, 

/ £dx- £dx= a= d{pdt) = - d(tdp) 

= - [ t(p)dp(p)+ [ (t(p)+T(p))dp= [ T(p)dp 

J a J a J a 



= T(£ )(£i - E ) + 0((E 1 - EoY), 

where T is the "rectangular domain" obtained by placing a "cut" a from ^{Eq) to 
7 and we get the (well-known) formula, 

^I(l(E))=T( 1 (E)). 

Since I(^(E)) only depends on E and not on the choice of "i{E), we get 2). 

Let p e be as in the introduction, and let q be defined in (1.16). Let G(x,£) be 
an analytic function defined in a neighborhood of J9 _1 (0), such that 

H p G = q-{q), (2.1) 

where we recall that (q) is the trajectory average, defined in (1.15). 

We will replace T*M by the new IR-manifold A eG = exp {ieH G ){T* M) (defined 
in a complex neighborhood of r ). Writing A e c 3 = exp (ieHa)(y,r)), we use 

p = (y, i]) as real symplectic coordinates on A e c- By Taylor expansion, we get 

p e (exp(iif G (p))) = ( P + ieq)(exp(ieH G (p)) + 0(e 2 ) = (2.2) 
p{p) + ie(q - H p G)(p) + 0{t 2 ) =p + ie{q) + 0(e 2 ). 

Recall the assumptions (1.17), (1.22), where we shall assume for simplicity that 
Fq — 0. (This is no real restriction, since we can always replace p e by p € — ieF$.) 
Since the Poisson bracket {p, Re(g)} is zero, we see that every component of the 
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set A 0) o = {p = 0, (q) = 0} is a smooth Lagrangian torus. Assume for simplicity 
(as in the introduction), that we only have one such component. Near this torus, 
p, Re (q) form an integrable system, so we can find a real and analytic canonical 
transformation hT x from a neighborhood of Ao,o to a neighborhood of £ = in 
T*T 2 , so that poK and Re (q) o k (and hence also (q) o k because of (1.17)) become 
functions of £ only. Here T 2 = (R/2ttZ) 2 . 

We can do this in the following way: Let k E ,F be the Lagrangian torus given by 
p = E,Re{q) = F, for (E, F) G neigh (0, R 2 ). ' Let 7i(£,F) be the cycle in A EiF 
corresponding to a closed ifp-trajectory with minimal period, and let 72 (E, F) be a 
second cycle so that 71,72 form a fundamental system of cycles on the torus A E>F . 
Necessarily 72 maps to the simple loop given by (q) = F in the abstract quotient 
manifold p~ 1 (E)/HH p . Now it is classical (see Arnold [Ar]) that we can find a 
real analytic canonical transformation k : neigh (r/ = 0,T*T 2 ) 3 (y,rj) 1— > (x, £) G 
neigh (A , , T*R 2 ), T 2 := (R/2ttZ) 2 such that 



Vj 



±([ idx- [ £dx), (2.3) 

27T J lj{E ,F) J ^(0,0) 



where E, F depend on (x, £) and are determined by (x, £) G A^, i.e. by = 

p(x,0,i 71 = Re(g)(a:,0- (See also [9].) 

In the following we sometimes write p instead of p o k and similarly for (q) (cf 
(1.17)): 

p = K0> (q) = (qK0- 

Then if p = J^af - - From (2.3) and the discussion at the beginning of this 
section, we see that p = p(£i) in the new coordinates, so 

# P = c(6)^-, P = p(6),c=J|^0. (2.4) 
The assumption (1.22) implies: 

dp dRejq) 

Thus 

p e =p(6)+^(g)(0+r e (x,0, (2.6) 

where 77 = C(e 2 ) and p, (q) satisfy (2.5). 

Now look for a " Lagrangian" torus T in the complex domain of the form 

£ = <f>\x\ x G T 2 , (2.7) 
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with <fi grad-periodic (in the sense that V0 is single-valued on T 2 ) and complex- 
valued, (j)' = 0(e), £<£«1, such that p e \ = 0. We get the eiconal equation 

P(^) + ie{q) (<&) + r e (x, cf>' x ) = 0, (2.8) 

where r t = 0(e 2 ). Write pfa) = <i + O^i), <?)(0 = <i + + 0(£ 2 ), c G R, 
c, Re 6 7^ so that: 

((c + ^e)— + ^>^r)0 + ^(s, <&) = 0, 

where 

F e (x,o = 0( e 2 + ^ 2 +ei). 

For notational convenience, assume that a = 0, b, c = 1. Look for = e0, with 
i>' = (9(1), e < e < 1. Then we get 

+ k— )V + G e , ? 0r,O = 0, ( 2 - 9 ) 



(9X1 5^2 

with 

Let H m (T 2 ) denote the standard Sobolev space of order m. In the following 
estimates, m > 1 is fixed. Using a standard result about non-linear functions of 
Sobolev class functions (see [1]), we get 

1) If (e-^AjV = 0(1) in H m , then G^(x,ip' x ) = 0{e{\ + e)) in # m . 

2) If {e-^dMj = Oil) in # m for j = 0, 1, then, 

|| [G^x, mi\\ H ™ = 0(e(i + e))\\(-d Xl , d X2 ) - ^o) Ik- 



3) If 



then 



(— h ?e— — )w = v, with -u, f periodic, 

OX i OX2 



C 

(e" 1 a ri ,a r2 )u||H'» < — Iklk" 



We shall find solutions to (2.9) that are grad-periodic functions of the form 

ip = "4> P er + a(exi + ix 2 ) + b(exi - ix 2 ), ( 2 -10) 
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with a given complex constant a = 0(1), and where the periodic function ip per and 
the complex constant b will depend on a. If Fuik) denotes the Fourier coefficient 
of u at k, we get the system: 

2eb + F(G e t(x, ip' per + a{ex l + ix 2 )' + b{ex 1 - ix 2 )'))(0) = 

d d , 

(- h ie- — )-^per + G €> e{x, ip pcr + a(exi + ix 2 )' + &( ea; i - ix 2 )') + 2eb = 0. 

C/Xx UX 2 

(2.H) 

We will find the solution as a limit of a sequence 

^(j) = ^g) + a ( eXi + ia . 2 ) + ftO)( eXl _ ^ 2 ) ; j = o, 1, 2, ... 

with V (0) = a(e^i + ix 2 ), (and &(°) = 0, rf$l = 0) where we impose 

(^ + + GetM^') = 0- 

The last equation gives the following system analogous to (2.11) that we label (Sj): 

2e&0-+i) + JF(G ei? (x, ^2' + a(en + ix 2 )' + b^(e Xl - ix 2 )'))(0) = 

" ! ' ^^-)^ 1) + G e f(x, ^<£' + a(ex! + ix 2 )' + 6^(^i - ix 2 )') + 2S 1+l) = 0. 



v 9xi &V per ' pcr 

From (S ), and the facts 1), 3), we get 



|6«| =<9(l)(|+c), 



ii(V,^ffiiifl-<o(i)(i+g. 

e e 
For j > 1, we consider (Sj) — (Sj_i) and get, using also 2), 

| 6 (j+D _ 6 (»)| < e)(||( V,0xjtyg> - l&^llfl- + \b U) ~ 



- e d Xl ,d X2 )(^-^ r )\\H^< 



per rper 

c c 

This implies that 



| 6 0-+D _ 6 0)| + ||(-^,^ 2 )(^ - Vg)ll^ < (0(l)U + e))'" +1 , 



per ' !»•!•' ; 
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and since e <C e <C 1, we see that the schema converges towards a solution to (2.9) 
of the form (2.10), with a = 0(1) (given), and 

\b\ + 11(^,^)^11^ = 0(l)(|+e). (2.12) 

For we then have 

= P er + ea(exi + ix 2 ) + e&(erci — ix 2 ), (2-13) 

with ea = 0(e) given, \eb\ + ||(e -1 9 xi , d X2 )(f) per \\ H m = 0(l)(e + e 2 ). In particular, 

= aee + C(e(e + e 2 )), |^ = ia? + C(e + e 2 ). (2.14) 



In this discussion m is fixed and the estimates are uniform with respect to e. 
Clearly only depends on the choice of ea (with m being fixed). As in [14], we see 
that depends holomorphically on ea, and extends holomorphically in x to some 
(e, e)-dependent domain in such a way that the dependence of ea is still holomorphic. 
In the preceding constructions, everything works the same way, if we replace T 2 by 
T 2 + iy, \y\ < 1/C, so it follows that <fi extends in x to a complex neighborhood of 
the real torus, which is independent of e, e, and that the preceding estimates remain 
valid here. 

Write = a , when e is fixed. Let Tf £ = 0'(x), x G T 2 . Let ij(r^), j = 1,2 
be the corresponding actions with respect to £idxi +£ 2 dx 2 . From (2.12), (2.13) (or 
simply (2.14)) we get: 

7i(r ) = 27r?e(a + 6) = 27r?e(a + o(L + e)), (2.15) 
J 2 (r^) =2irie{a-b) = 2ne(ia + + e)), 

We are interested in finding a such that both actions are real. This leads to 

Im (a + 0(1 + e)) = 0, Im (ia + C(| + e)) = 0, 



(2.16) 



i.e. 

Rea + 0(i + 7) =0 
lma + 0(| + e) = 0. 

Here the O-terms are real parts of holomorphic functions, so they remain 0(e + e/e) 
after derivation with respect to Re a, Ima. By the implicit function theorem, we 



11 



therefore have a unique solution to (2.16), which is 0(e + e/e), and correspondingly 
ea = 0(e + e 2 ). Recall that ea is independent of the choice of?, so if we take e — y/e, 
we get 

ea = 0{e). (2.17) 

For this particular <f>, we have 

d Xl( j> = 0(e 2 ), d X2 (j> = O(e) in H m . (2.18) 

If we do not make the simplifying assumption that J|j-(0) = 1, ||^(0) = 0, 

^■(0) = 1, then the earlier discussion tells us that we have solutions of the type 
(f) — eip, — ip per (x) + aa(x) + b/3(x), with 

«W = e^^(0)xi + i — (0)x 2 , 

of \ d (<l) / n x .9{p + ie(q)) f , 

Observe that if we put, 



z _ = djp + iejq)) ^ d + diejq) ^ d 



d£x dx 1 <9£ 2 9x2 

then 

Za = a , z , 3 = 2f ^±ilM> (0) f) (o) + o. 

The earlier discussion goes through without any changes. Especially, in the case 
a = 0, the corresponding is independent of e. 

Let now ( vary in neigh (0, C 2 ). Put z(() = p((i) +^ e (l)(C)- Then the discussion 
above can be applied with p e (x, £) replaced by 

p e (x, c + - z(0 = p(Ci + 6) - p(Ci) + ie«?XC + - (?>(0) + <^ 2 )- 
We get a solution to the eiconal equation 

Pe (x,c + <)-2(C) = o 

of the form 

^(x,c) = ^pe r (*,C) + HCM*,C), 

where 
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depending holomorphically on (. (So we choose a = in the earlier discussion, 
but compensate for this by introducing a ^-dependence and even varying the energy 
level z(C).) 

As before, we get 

\\(- e d Xl ,d X2 )ij pei \\ H m = 0(e), \b\ = (9(e), (2.19) 

uniformly with respect to (. Moreover, since the problem depends holomorphically 
on (, it is easy to see (for instance by working in a space of holomorphic functions 
of ( with values in H m ) that V^, b depend holomorphically on (. Notice that 

4>(xX)--=x-( + ^(xX) (2.20) 

solves the eiconal equation 

Pe (x,dJ(x,C))-z(O=0. (2.21) 

Write b(()(3(x, () + x ■ ( = x ■ 77, where 77(C) depends holomorphically on ( and 
satisfies 

Vl (()=( 1 + 0(6 2 ), m (0=C2 + O(e). 

Let ((■!]) with Ci(v) — Vi + ^( e2 )> Giv) = V2 + 0(e) denote the inverse. Then with 

(f) pe r(x : rj) = VV-O^C), we have 

$(x, C) = x ■ rj + (f) pcr (x, rj) =: (j)(x, rj), (2.22) 

solving 

p e (x,d x <l>(x,ri))-p e (ri) = 0, p e (r)) = z(((r))) = pfa) + ie{q)(rj) + C(e 2 ), (2.23) 
while (2.19) gives 

\(^d xl ,d xa )<f> vet (x,ri)\ = 0(e), (2.24) 

for 1 in a fixed complex neighborhood of T 2 and as usual, we get corresponding 
estimates for d"d^(p per (x,r]) from the Cauchy inequalities. We normalize the choice 
of 4> per (x, i]) by requiring that 



(0per(-, V)) = J^y 0per(z, T])dx = 0. 



Then 

«e : v),ri)*-> (x, (j)' x (x, r])) (2.25) 
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maps a complex (e-independent) neighborhood of the zero section of T*T 2 onto 
another neighborhood of the same type (containing an e-independent neighborhood 
of f = 0). (2.23) shows that 

p e o K e = p e . (2.26) 

By construction, we also know that K t conserves actions along closed curves. 

Using that (f)' v (x,r]) = x + 0(e), (f>' x (x,r)) = i] + 0(e 2 ,e) together with (2.24), 
which also holds with per replaced by its gradient, we see that 

« e (y, V) = (V + O(e); m + 0(e 2 ), m + 0(e)). (2.27) 
In particular, we have 

Imx = (9(e), Imfi = (9(e 2 ), Im£ 2 = (9(e), (2.28) 
on the image of T*T 2 . We can therefore represent K e (T*T 2 ) by 

Imx = G^(Re(x,0), Im£ = -G^(Re (x, £)), (2.29) 

where G is a smooth, a priori grad-periodic function which satisfies, 

d ( G, d X2 G = (9(e), d Xl G = (9(e 2 ). (2.30) 

Since k € conserves actions, the actions along closed cycles in K e (T*T 2 ) are real and 
it follows that G is single- valued. We may assume that G = (9(e). Let x(0 be a 
standard cutoff around £ = and let M e be given by 

Imx = G'z(Re(x,0), Im£ = -G^(Re (x, 0), (2.31) 
where G(Re(x,0) = x(ReOG(Re (x, £)). 

Then M e is an IR-manifold which coincides with T*T 2 outside a (complex e-indepen- 
dent) neighborhood of £ = 0. Moreover, we know that M t is an e-perturbation of 
T*T 2 , along which we have 

Im& = - X (Re0^ 1 (Re(x,0) = <?(e 2 ). 

It follows that outside the neighborhood of £ = 0, where M t coincides with K e (T*T 2 ), 
we have 

\tePe^\ + \\*n P ett)>± (2.32) 

Now recall the initial global situation, that we simplified the original principal 
symbol by composing with expieHc for the function G in (2.1) and then further by 
k, introduced prior to (2.4). 
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We introduce an IR-deformation M e of real phase space which is anjs-deformation, 
equal to real phase space away from r , and equal to expieH G o k(M € ) near T = 
p _1 (0) fl T*M. Then we have achieved the following: 

Proposition 2.1 a) There exists an analytic real canonical transformation K e : 
neigh (f = 0, T*T 2 ) -> neigh (exp (ie# G )(A , ), M e ), such that 

p e o K e =p £ (rj), (2.33) 

where p e is given in (2.23). 

b) Away from the small neighborhood, where (2.33) holds, we have 

l^, M J + i|Imp e|M J>i (2.34) 

Here p e denotes the original principal symbol of the perturbed operator. 

It is now clear that the main result of [MeSj] can be applied to give the full 
asymptotics for all eigenvalues of P t in a domain |Rez| < 1/(9(1), \lmz\ < e/0(l), 
for e > small enough and for h < h(e) > small enough depending on e. It is 
not apriori clear however what kind of uniformity with respect to e we may have 
in this result. We shall employ quantum Birkhoff normal forms in the next section 
and obtain a more uniform result, valid for e > h s for any fixed S e]0, |[. 

3 Formal spectral asymptotics 

As in [9] (see also [14]) we can implement k € by an elliptic Fourier integral operator 
U — U e : L|(T 2 ) — > H(M e ) which is microlocally defined from a neighborhood of 
f = in T*T 2 to a neighborhood of exp ieH G (A 0fi ) in M e . Here S = (S u S 2 ) E R 2 , 
with Sj = f^.^dx, and jj = 7^(0,0) are introduced prior to (2.3). Ls(T 2 ) denotes 

the space of locally L 2 -functions u on R 2 satisfying the Floquet periodicity condition: 

u(x - 7 ) = e&<* s+ % a °\ 7 e (2vrZ) 2 , (3.1) 

where a = (a®, a®) G Z 2 is a Maslov index. By abuse of notation, we still denote 
by P e , the conjugated operator U~ 1 P e U e . 

We have an analytic /i-pseudodifferential operator P e on T 2 (defined microlocally 
near £ = 0), of order in h, with leading symbol independent of x: 

Pe(0=P(ti)+K(q)(0 + O(e 2 ), (3.2) 
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defined in a fixed complex neighborhood of £ = in T*T 2 , depending holomorphi- 
cally on e G -D(0, e ). The full symbol is 

oo 

P e (x,Z;h)=Y t h i Pj (x,Z,e), (3.3) 
j=o 

with £, e) holomorphic with respect to (x, £) in a j-independent complex neigh- 
borhood of £ = and C°° with respect to e G [0, e [, with p (x,£,e) = p t {£)- 
Following the standard Birkhoff normal form procedure, we shall remove the x- 
dependence in the pj by means of conjugation by an elliptic /i-pseudodifferential 
operator of order 0. Let A be an /i-pseudodifferential operator of order 0. Recall 
that 

e A Pe- A = e adA P = ^ ad A R 
Let the full symbol of A be of the form YlT=o h k au- Then on the operator level, 



oo oo oo 



£=0 fc=0ji=0 j k =0 

oo 

?1=0 

with s = p , si = \H ao p +pi = iH po a + p 1 , s n+l = iH Po a n + s n+1 ,..., where 
s n+ i only depends on a , a n _i and is the sum of the coefficients for h n+1 from the 
terms 

^ Jl+ -^ +m (^ad an )...(iad a ,J(p,), 

with 

ji + ••• + jk + ? + k < n + 1, ji, j fe < n, or k = 0, £ = n + 1. 
Notice that 

and that we can solve 

H po a = b(x,0-(K;0), ^T 2 , (3.4) 

with II < 0(l)e 1 ||6||//m. As already noticed in the preceding section, the 

same equation can be solved in a complex domain {x G T 2 ; |Imx| < C2}, and we 
get 

sup \a(x,Z)\< C{Cl ' C2) sup \b(x,0\, (3.5) 

\lmx\<C2 e \lmx\<C\ 



16 



if C\ < C 2 . The shrinking of the domains in (3.5) is not a problem, since we can 
take a sequence of such domains with Cj \ > 0. 

By solving equations of the type (3.4), we can determine do, Oi, ••• successively, so 
that Sj = Sj(C,, e) are independent of x. Assume by induction that Va, = G(e~ 1 ~ 2j ), 
for j < n — 1 (in a complex domain, so that we have the same estimates on the 
derivatives of Vaj). Then the general term in s n+1 is 

0(l)e~ 1 ~ 2j ' 1 ...e~ 1-2j ' fc = 0(l)(-) 2( - jl+ - +jk)+k . 

Here, 

2(ji + •• + 3k) + k = 2(j\ + .. + j k + k) - k < 2(n + 1 - £) - k = 2n + 2 - 21 - fc. 

So this quantity is (9(l)(-) 2ri except possibly when 2£ + k < 2, i.e. when fc = £ = 
or when k — 1, I — 0. In the first case we get the coefficient for in p which 
is 0. In the second case, we get the coefficient for h n+1 in h^ 1+1 (^ad aj ){po) with 
j\ < n, which is 0(l){\) 1+2 i\ Here 1 + 2ji < 2n. Thus s n+ i = C( e - 2n ) (in a 
complex domain). We can choose a n periodic, with iH po a n = — s n+1 + (s n +i( - ,£)) 
and with Va n = (9(e _1_2n ). This completes the induction step and we find a k with 
Vdfc = 0(e~ l ~ 2k ) in a fixed complex neighborhood of T 2 x {£ = 0} such that if 

N-l 

AW = J2 h k a k , 

k=0 

then 

oo 
n=0 

where Pn^(^,e) = 0{e~ 2< " n ~ r>+ ) and = pi 00 -* is independent of x and N, for 
n < N. From this we get the following formal spectral result: 

Theorem 3.1 Under the assumptions above, there exists a constant C > such 
that if 5 > is fixed and h^~ 5 < e < 1/C, and < h < h(5) with h(5) > small 
enough, then in the region 

_ 1 llmzl 1 . 

\Rez\<-, L T A <£, (3-7) 

P has the following quasi- eigenvalues : 

oo „ o 

**~£ ^(M* - ^ - y)- e), k e Z 2 . (3.8) 

ra=0 
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Here, S G R 2 , a G Z 2 were introduced in the beginning of this section, andp^ = p e 
is given in (3.2). 

We leave undefined, the notion of quasi-eigenvalue, and interpret the above theo- 
rem as the formal consequence of the reductions above and the fact that the functions 

e k (x) = e ix < k -^~^\ keZ 2 , 
form an orthonormal basis in L|(T 2 ). 

4 Justification via a global Grushin problem. 

In this section we outline how Theorem 3.1 actually gives all eigenvalues in the 
rectangle (3.7). As in [14], [9] we construct an auxiliary, so called Grushin problem. 
Actually, this construction is identical with the one in [9], so we shall only recall the 
main steps. 

For C > sufficiently large, let I(C, e) (depending also on h) be the set of all 
k G Z 2 , for which the values z k in (3.8) belong to the rectangle (3.7). Recall that z k 
correspond to the orthonormal family of functions e^, defined after Theorem 3.1. 

Let n e ,M e be as in Proposition, 2.1 and let U e be the Fourier integral operator 
quantization of n e introduced in the beginning of Section 3. With defined there, 
let A be a natural asymptotic limit. Define 

R+ : H(M e ) -> C I(c ' e \ (4.1) 

by 

R+u(k) = (e^- 1 ^)^. (4.2) 

Notice that R + is a globally welldefined operator modulo some indetermination of 
norm 0(h°°), since e A U~ 1 u is microlocally welldefined in a neighborhood of the zero 
section in T*T 2 . Similarly, we define R- : C I(c ' e) -> H(M e ), by 

R.u^ = u4k)U e e- A e k . (4.3) 

keI(C,e) 

Then for z in the rectangle (3.7), with an increased value of C, the problem 

(P - z)u + P_m_ = v, R + u = v+, (4.4) 

has a unique solution (u,U-) G H(M e ) x Q 1 ^) for every (v,v+) G H(M e ) x C I<yC ' e \ 
(Here we assume for simplicity that P is a bounded operator, otherwise we would 
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have to work with modifications of H(M e ) of Sobolev type, depending on additional 
order functions. See the appendix in [9] for more details and further references.) We 
have the corresponding apriori estmate 

C 

IMI + IMI < — (H +e|K||), (4.5) 
and if we write the solution 

u \ ( E E+ \ { v \ (4 6) 



then modulo 0(h°°), E_ + is the diagonal matrix ((z — Zk)5j t k), where Zk are given 
in (3.8). 

Recall from [9] that the verification of these facts consists of half-estimates away 
from A 0j o and the exploitation near A 0j o of the reduction to a translation invariant 
operator on T 2 in the preceding section. Since the eigenvalues of P in our rectangle 
are precisely the values z for which E_ + (z) is non invertible, we get 

Theorem 4.1 Under the assumptions of Theorem 3.1, there exists a constant C > 
such that if 5 > is fixed and h^~ 5 < e < 1/C , and < h < h(5) with h(5) > 
small enough, then in the region 

|R.z| <-,!_!<-, (4.7) 
the eigenvalues of P are simple and given by 

**~£ h»$r\h(k - — - -), e),k G Z 2 , (4.8) 

ra=0 



with one eigenvalue for each k such that z^ belongs to (4-7)- Here, S G R 2 ; a G Z 2 
were introduced in the beginning of Section 3, and the p£°^ were constructed prior 
to Theorem 3.1. Further, pfi°\t, e) = p(&) + ie(q)(£) + 0(e 2 ). 



5 Application to barrier top resonances. 

We extend the domain of validity of one of the results of section 7 in [9], by using 
Theorem 4.1 as the new ingredient. The discussion that follows will therefore only 
be a brief recollection of a part of Section 7 in [9], and we refer to that work for 
more details. Let 

P = -h 2 A + V(x), p(x,0 = e + V(x), (x,OeT*R 2 = R\ (5.1) 
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satisfy the general conditions for defining resonances near the energy level E > 0. 
Assume that V(0) = 0, W(0) = 0, V"(0) < and that V is everywhere analytic. 
After a linear change of x-coordinates, we have near x = 0: 

2 A 

p(x,0 -E = Y1 - x )) +P3(x) +p 4 (x) + (5.2) 

i 

where \j > and p v is a homogeneous polynomial of degree v. Also assume that 
(0,0) is the only trapped point for the Hp-How on the real energy surface p~ 1 (Eq). 
We assume A = (Ai, A 2 ) fulfills the resonance condition 

A • k = 0, for some O^keZ 2 . (5.3) 

Somewhat roughly, the problem of determining the resonances near Eq is then equiv- 
alent to determining the eigenvalues of P — Eq near 0, after the change of variables, 
x = e t7T ^x (and £ = e _M "/ 4 £) near 0, and we get a new operator with symbol 

-i(p 2 (s, + ie 37ri/4 p 3 (£) + ie^ i/A p±{x) + ...) = £), (5.4) 

i 

Dropping the tildes for the new variables, we are then interested in eigenvalues E 
of Q = g(ar, /lAr) with |£| ~ e 2 , /i" 5 < e < 1, < 5 < 1/2. Write x = ey, h = h/e 2 . 
Then hD x = ehD y and 

e- 2 q(x, = e- 2 g(e(y, 77)) = p 2 (y, 77) + 2ee 3 ^ 4 p 3 (l/) + ^(e 2 ), 

in a region | (y, 77) ] = 0(1), where the corresponding eigenfunctions are concentrated. 

The resonance condition (5.3) implies that the H P2 -fiow is periodic with a period 
T > 0, independent of the energy level. Using Theorem 4.1 in the discussion of 
section 7 in [9], we get the following variant of Proposition 7.1 of that paper: 

Proposition 5.1 Let (p 3 ) denote the average of P3 along the trajectories of the 
Hamilton vector field of P2 in (5.4), and assume that (p 3 ) is not identically zero. Let 
F e R be a regular value of cos(37r/4)(p 3 ) restricted to p 2 " 1 (l), and assume that T 
is the minimal period of the H ' ^-trajectories in the torus A 1>Fo given by 

Ai,fb : P2 = 1, cos (£) (p 3 ) = F . 
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Let e satisfy 

Then for z in the set 
1 



h s « e < e , < e < 1, < 5 < -. 



(5.5) 



1 1 
1 + 



0(1)' 0(1). 



+ ie 



0(1 



0(1) 



, k e z 2 



the resonances of the form E — ie 2 z are given by 

z = P^(k-j)-^,e,h^+0(h 00 ), h= t _ 

fmt/i precisely one resonance for every k). Here P , e; /i) has an expansion as 

CO 



n=0 



where 



The coordinates £x = £i(E) and £ 2 = £,2{E,F) are the normalized actions of 

Ae,f ■ P2 = E, cos(^) (p 3 ) = F, 



for E e neigh(l,R) ; F e neigh(F , R) ; given by 



1 

2^ 



77 - 



7i (£,F) 



vdy }, j = 1, 2, 



(5.6) 



7,-(i.Jb) 



to£/j jj(E,F) being fundamental cycles in Ae,f, such that ji(E,F) corresponds to 
a closed H ' ^-trajectory of minimal period T. Furthermore, 



= vdy, j = 1,2, S = (Si, £ 2 ), 



(5.7) 



and a E Z 2 is fixed. 

The interest of this result (as well as of Theorem 4.1) compared to the corre- 
sponding ones in [9] is that we can reach small but /i-independent values of e. On 
the other hand our method does not immediately seem to be able to hand as small 
values of e as in [9], and the results there give a desrciption of how the negative 
powers of e appear in our estimates of the terms in the asymptotic expansion of the 
symbol P(f,e; h). 



21 



References 



[1] S. Alinhac, P. Gerard, Operateurs pseudo-differentiels et theoreme de Nash- 
Moser, Savoirs actuels, InterEditions/Editions du CNRS, Paris 1991. 

[2] M. Asch, G. Lebeau, The spectrum of the damped wave equation, preprint, 1999. 
See http://www.math.u-psud.fr/~biblio/rt/1999/. 

[3] D. Bambusi, S. Graffi and T. Paul, Normal forms and quantization formulae, 
Comm. Math. Phys. 207 (1999), 173-195. 

[4] L. Boutet de Monvel, V. Guillemin, The spectral theory of Toeplitz operators, 
Annals of Math. Studies 99, Princeton University Press, 1981. 

[5] Y. Colin de Verdiere, Sur le spectre des operateurs elliptiques a bicar- 
acteristiques toutes periodiques, Comment Math. Helv. 54 (1979), 508-522. 

[6] B. Helffer, D. Robert, Puits de potentiel generalises et asymptotique semi- 
classique, Ann. Inst. H. Poincare 41 (1984), 291-331. 

[7] B. Helffer, J. Sjostrand, Resonances en limite semiclassique, Mem. Soc. Math. 
France (N.S.) 24-25 (1986). 

[8] M. Hitrik, Eig en frequencies for damped wave equations on Zoll manifolds, 
Asymptot. Analysis 31 (2002), 265-277. 

[9] M. Hitrik, J. Sjostrand, N on-self adjoint perturbations of selfadjoint operators 
in 2 dimensions I, Preprint 2003. 

[10] M. Ikawa, On the scattering matrix for two convex obstacles, Hyperbolic equa- 
tions and related topics (Katata/Kyoto, 1984), 63-84, Academic Press, Boston, 
MA, 1986. 

[11] V. Ivrii, Microlocal analysis and precise spectral asymptotics, Springer- Verlag, 
Berlin, 1998. 

[12] N. Kaidi, P. Kerdelhue, Forme normale de Birkhoff et resonances, Asymptot. 
Analysis 23 (2000), 1-21. 

[13] G. Lebeau, Equation des ondes amorties, in Algebraic and Geometric Methods 
of Mathematical Physics (Kaciveli 1993), 73-109, Math. Phys. Stud., 19 Kluwer 
Acad. Publ., Dordrecht, 1996. 



22 



[14] A. Melin, J. Sjostrand, Bohr-Sommerfeld quantization condition for non-self- 
adjoint operators in dimension 2, Asterisque, to appear. 

[15] G. Popov, Invariant tori, effective stability, and quasimodes with exponen- 
tially small error terms. II. Quantum Birkhoff normal forms, Ann. Inst. Henri 
Poincare, Phys. Th. 1 (2000), 249-279. 

[16] J. Sjostrand, Resonances associated to a closed hyperbolic trajectory in dimen- 
sion 2, preprint, Autumn 2002. 

[17] J. Sjostrand, Asymptotic distribution of of eig en frequencies for damped wave 
equations, Publ. Res. Inst. Math. Sci. 36 (2000), 573-611. 

[18] J. Sjostrand, Semi-excited states in non- degenerate potential wells, Asymptot. 
Analysis 6 (1992), 29-43. 

[19] J. Sjostrand, Singularites analytiques microlocales, Asterisque 85 (1982). 

[20] J. Sjostrand and M. Zworski, Asymptotic distribution of resonances for convex 
obstacles, Acta Math. 183 (1999), 191-253. 

[21] A. Weinstein, Asymptotics of eigenvalue clusters for the Laplacian plus a po- 
tential, Duke Math. J. 44 (1977), 883-892. 



23 



